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Abstract 

DNA molecule conformation dynamics and stretching were made on semi-circular surfaces with different radii (500 
to 5,000 |jm) in microchannels measuring 200 |jm x 200 |jm in cross section. Five different buffer solutions - Ix 
Tris-acetate-EDTA CTAE), 1 x Tris-borate-EDTA (TBE), Ix Tris-EDTA (YE), Ix Tris-phosphate-EDTA (TPE), and Ix 
Tris-buffered saline (TBS) solutions - were used with a variety of viscosity such as 40, 60, and 80 cP, with resultant 
10""* < Re < 10"^ and the corresponding 5 < Wi < 12. The test fluids were seeded with JOJO-1 tracer particles for 
flow visualization and driven through the test channels via a piezoelectric (PZT) micropump. Micro particle image 
velocimetry (|jPIV) measuring technique was applied for the centered-plane velocity distribution measurements. It is 
found that the radius effect on the stretch ratio of DNA dependence is significant. The stretch ratio becomes larger 
as the radius becomes small due to the larger centrifugal force. Consequently, the maximum stretch was found at 
the center of the channel with a radius of 500 |jm. 
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Background 

Fully stretched DNA molecules are very important with 
regard to advancing the genomic sciences and analyses 
in order to understand the physical and biological proper- 
ties of DNA, including the ability to directly manipulate 
and visualize single DNA molecules. In fact, engineering 
DNA stretching would be a key step in the develop- 
ment of the next generation of biological microfluidic 
devices [1,2]. 

Microfluidics is the study of behavior manipulation 
and control of fluids confined to micrometer dimensions, 
typically 1 to 100 \im. Transport in the microchannels is 
the major phenomenon; it includes flow detections, liquid 
transport, control of molecular transport lil<e DNA mol- 
ecule conformation dynamics, measurement of buLk-level 
rheological properties, and separation techniques with bio- 
physical and genomic applications because they generate 
defined fluid flows that manipulate large DNA molecules 
[3]. In addition, understanding the complex behavior of 
DNA molecules flowing in microchannels is essential 
to the realization of lap-on-a-chip (LOG) and micro 
total analysis system (j^TAS) intended to systematically 
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manipulate, process, and analyze these molecules. The 
presence of DNA molecules gives the fluid viscoelastic 
behavior that may change the base flow pattern in 
curved channels [4] . 

Two general approaches to DNA stretching are in 
common use: DNA is stretched in a solution as it flows 
through a microchannel, or it is stretched on a solid 
surface. For the latter, the conditions required for sig- 
nificant DNA stretching include high shear rates and 
high pressure gradient operations with a pressure-driven 
flow, due to non-slip boundary conditions on the wall. 
The shear flow existing at the channel walls could stretch 
DNA molecules. The degree of stretching is correlated 
with the Weissenberg number of the flow, Wi = r where t 
is a characteristic relaxation time for the molecule in the 
solution and is a characteristic shear rate based on the 
flow in the channel. For the past two decades, DNA mole- 
cules have served as a model system for single molecule 
semi-flexible polymer (larger persistence length of ap- 
proximately 65 mm) dynamics and can be fluorescently 
labeled for direct observation with videomicroscopy [5], 
revealing a DNA solution non-equOibrium microstructure, 
DNA- solvent interactions, and DNA macromolecular 
transport phenomena. Moreover, increased interactions 
between DNA molecules and channel surfaces result in 
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non-Newtonian flow behavior, even in a dilute DNA 
molecule solution. 

When the laminar flow passes through the curved 
channels/ducts, the centrifugal force pushes the fluid 
from the center of the channel when the bulk fluid flows 
with high velocity to the outer side, while the fluid at the 
outer wall is pressed either upwards or downwards, thus 
producing two vortices to fill the entire channel at a 
cross section along the downstream [6]. The mean flow 
velocity and the curvature of the channel can determine 
the centrifugal force, which is governed by an important 
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Figure 1 Fabrication of the present curved channel. 



dimensionless parameter of Dean number (Dn = Re (dh/ 
R)'^'^), including the flow Reynolds number (Re) and the 
duct hydraulic diameter ((i^) to the curvature of the 
channel/duct (R). Here, the Reynolds number is defined 
as Re = pUdh/fi{y), where p is the solution density, U 
is the average velocity, ^ is the solution viscosity, and 
(/) is the solution shear rate. 

Research on shear flow [3,7,8] has been conducted in 
order to model the conformation of DNA molecules 
for an extended time. These studies reported that the 
stretching of the DNA molecules subject to shear flow 
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Table 1 Fabrication parameters of curved microchannel 


Relevant parameter 


SU-8 fabrication processes 


Spin coating 


550 rpm (80 s) 




800 rpnn (70 s) 


Soft bal<e 


65°C (3 min) 




95°C (21 min) 




65°C (3 min) 




Room temperature (30 min) 


Exposure 


Total time 30 s 


Post exposure bal<e 


65°C (3 min) 




95°C (12 min) 




65°C (3 min) 


Hard bal<e 


95°C (3 min) 


Remarl<s 


Substrate type: silicon wafer 




Photoresist: SU8-2100 (MicroChem) 




Depth: 200 |jm 




Photomask: film mask 




(FUJI HPB-S 7mil, 20,000 DPI) 


PDMS fabrication processes 


PDMS prepared 


10:1 Sylgard-184 A/Sylgard-184 B mixture 


Bal<e 


70°C (21 min) 



is a function of Wi and r in the flow. In this study, 
XDNA molecules were stretched on curved wall sur- 
faces in different curved ducts in pressure-driven flows 
and visualized as well as measured via micro particle 
image velocimetry (|iPIV) and an optical system. Special 
attention will be paid to examining the effect of different 



radii of the curved duct (i.e., Dn), buffer solutions, and 
the viscosity of the solution. Moreover, viscoelastic (i.e., 
non-Newtonian) flow in dilute DNA solution will also 
be examined. 

Methods 

PDMS flow cell fabrication 

In this study, we used a 200 i^m x 200 |im microchannel, 
as shown in Figure 1. The polydimethylsiloxane (PDMS) 
channels were fabricated in-house at the University 
Microsystem Laboratory by casting open ten concentric 
circular-slot channels from PDMS and sealing it with 
the same material. At the center of these ten concentric 
circles, an up/down plenum was drilled to allow the buf- 
fer solution and DNA molecules to flow in/out; thus, the 
circular ducts became two symmetric half semi-circle 
ducts. The casting mold was made by SU-8 deep UV 
lithography. The detailed SU-8 mold design and PDMS 
curved channel fabrication created through UV lithog- 
raphy can be found in [1], with slight modifications for 
the photomask. Table 1 lists the fabrication parameters 
of the present curved microchannels. 

For the tested channels, precise information on their 
dimensions is extremely important to obtain an accurate 
evaluation of this microchannel. The depth, width, and 
length were measured optically within an accuracy of 
±0.2%. The surface roughness of the channel was mea- 
sured with a surface profilometer. During the experi- 
ments, the surface of the flow channel was so designed 
that the surface was kept hydrophilic in order to have 
the buffer solution flow through the microchannels with 
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Figure 2 Schematic for the present measuring instruments. 
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Table 2 Buffer solution used In the study 







1x TE 


1x TAE 


1x TBE 


1x TPE 


IxTBS 


Viscosity (cP) 


40 


60 80 


40 60 80 


40 60 80 


40 60 80 


40 60 80 


Sucrose (g/ml) 


1437 


1 .606 1 .726 


1 .437 1 .606 1 .726 


1.437 1.606 1.726 


1.437 1.606 1.726 


1.437 1.606 1.726 


Tris base concentration (mM) 




10 


40 


90 


90 


50 


EDTA concentration (mM) 




1 


1 


2 


2 


None 


Other ion concentration 


5.2 mM of liydrochloric 
acid 


20 mM of acetic 
acid 


90 mM of boric 
acid 


26 mM of phosplioric 
acid 


150 mM of sodium 
chloride 


pH 8 8 8 8 8 


Lambda DNA (ng/ml) 








0.0325 







JOJO-1 concentration (mM) 0.02 



a definite surface resistance. Pressure gradients in the 
present curved channels generated modified (due to cen- 
trifiagal force) parabolic flow, such that shear flow oc- 
curred near the channel walls. Furthermore, microfluidic 
semi-circular curved ducts created a periodic oscillating 
flow, in which flow pressure gradient alternated direc- 
tions at a definite time and extended observations of 
DNA molecules. 

DNA visualization and buffer solution preparation 

An experimental setup scheme combined with a laser 
light source (Ar-ion laser 488 nm/HeNe laser 532 nm) 
and scanning system used to implement |iPIV measure- 
ment is shown in Figure 2. The flow cell was mounted 
onto an epifluorescent microscope (IX71/FV300, Olym- 
pus, Tokyo, Japan) equipped with a x 40 magnification 
and NA 0.85 air immersion objective lens, following the 
description in [2,9]. The use of the i^PIV technique is 
very attractive in microfluidics because it helps to deter- 
mine the detailed flow phenomena of microsystems by 
utilizing flow-tracing particles to map the flow in the 
microchannels. Streak images and video microscopy as- 
sist in the investigation into the flow kinematics in the 
circular curved microchannels; jiPIV is used to quantify 
the flow field in the vicinity of the curved channels. In this 
study, the stained DNA molecules (JOJO-1, Invitrogen, 



Carlsbad, CA, USA) were used as seeding. The probe used 
to visualize the DNA was JOJO-1 at a dye with base pair 
ratio of 1:5. Incubation for the DNA and probe was initi- 
ated. The dyed A.DNA had a contour length (Ig) of 21 ^m 
and the longest relaxation time (r) of 7.6 s. 

Shear flow system 

A custom-made flow system was developed to enable 
the simultaneous generation of controlled shear flows 
and visualization of the DNA molecular conformation 
dynamics. The present DNA solution was found to be 
highly shear thinning at high shear rates, with a shear 
viscosity fi (cP) defined as the power law fi = 77/"" '' 
(one of typical relations). Flows of water and diluted 
DNA solution (ADNA, 31.5 x 10® D, 48,502 base pairs) 
are explored in the microchannels. The DNA solution 
was diluted to a concentration of 0.0325 pg/ml using five 
different buffers (see Tables 2 and 3 for details). In order 
to increase the viscous drag, the viscosity of the buffer 
solution was adjusted from 40 to 80 cP by adding a 
proper amount of sucrose. The test fluids, as stated pre- 
viously, were seeded with JOJO-1 tracer particles for 
flow visualization and driven through the circular curved 
ducts using a piezoelectric (PZT) micropump. A micro- 
filter was placed between the pressure regulator and the 



Table 3 Relevant parameters of the flow under study 



Parameter 




Value 




Pressure drop 




34 Pa, 44 Pa, 57 Pa 




Power consumption 




0.06 W, 0.068 W, 0.08 W 




DNA molecular concentration 




0.0325 pg/ml 




Working fluid viscosity, /j (cP) 


40 


60 


80 


Reynolds number. Re (xlO"^) 


1.2 to 1.87 


0.561 to 0.828 


0.326 to 0.486 


Dean number (xlO"'') 


1 .7 to 8.4 


0.8 to 4.1 


0.4 to 2.4 


Relaxation time, Tr (Rouse model) 


42 


631 


841 


Relaxation time, iz (Zimm model) 


3.1 


4.6 


6.1 


Relaxation time, T (present study) 


3.82 


5.6 


7.6 


Weissenberg number, Wi 


6.7 to 1 1 


7.2 to 1 1 .3 


8 to 12 



Hsieh ef al. Nanoscale Research Letters 2014, 9:382 
http://www.nanoscalereslett.eom/content/9/l/382 



Page 5 of 1 3 



350 



300 - 



250 



zi. 

^ 200 



^ 150 



S 100 



1 ' 1 ' 
■ ■ 1 X TE 


1 ' 1 ' 


1 


" • 1 X TAE 


|x = 40 cP ^ 








■ 1 X TBE 






- T 






^ 1 X TPE 






H = 60 cP 






__-f-^ — iT^locP 




I.I, 


I.I, 


1 



2.6 2.7 2.8 2.9 

Input voltage (V) 

Figure 3 Input voltage (DC) vs velocity for the present piezoelectric (PZT) micropump. 



3.0 



flow meter to eliminate any particles (>0.1 |im) or bub- 
bles (>0.1 \im). A tracing particle of stained DNA mole- 
cules was used for i^PIV measurements between the flow 
meter and the inlet and outlet of the channel. The mass 
flow rate was estimated through a stopwatch to count 
how long the buffer solution took to complete a flow loop, 



and the total weight of the buffer solution in a flow loop 
was measured by a microbalance. The mass flow rate 
found in this study was about 3 x 10" to 6 x 10" ml/min. 
The errors of the flow rate measurement were estimated 
to be less than ±3%. The DNA solution was delivered into 
the circular duct with two equal flow rate fluid delivery 
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Figure 5 DNA stretching and DNA molecule elongation, (a) Time history of DNA stretching at different Wi. (b) DNA molecule elongation 
length vs time. 



lines, with a very small Reynolds number in the range of 
0.326 X 10"^ to 1.87 x 10"^, in which molecular diffusion 
was a major mechanism for mixing. The Reynolds number 
was based on the shear rate-dependent viscosity ^, as 
stated previously. The characteristic shear rate y used 
for calculating Wi was taken to be the average velocity 
U divided by the channel half width w/2. 



[jPIV system 

The i^PIV utilizes flow-tracing particles (stained DNA 
molecules) to map the flow in the microchannels. The 
setup shown in Figure 2 was based on two-pulsed Nd: 
YAG lasers (New Wave SoloII, 30 mj, double cavity; 
New Wave Research, Inc., Fremont, CA, USA) firing 
on the second harmonic (green 532 nm). A detailed 
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Figure 6 Stretching ratio histogram for different buffers with different viscosities, (a) 40 cP, (b) 60 cP, and (c) 80 cP 
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Figure 7 Comparisons with the related previous studies for DNA stretching. 
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Table 4 DNA mean stretching rate 



Input voltage (DC) 


Buffer viscosity (cP) 


1x TE 


1x TAE 


1x TBE 


1x TPE 


1x TBS 


2.6 V 


40 


0.26 


0.252 


0.253 


0.265 


0.262 




60 


0.271 


0.266 


0.271 


0.2676 


0.2754 




80 


0.278 


0.283 


0.281 


0.28 


0.2844 


2.8 V 


40 


0.284 


0.2867 


0.283 


0.2867 


0.2922 




60 


0.288 


0.293 


0.289 


0.2917 


0.2953 




80 


0311 


0.301 


0.3 


0.3035 


0.308 


3.0 V 


40 


0.302 


0.309 


0.302 


0.3031 


0.3061 




60 


0317 


0.315 


0.307 


0316 


0315 



80 0.318 0317 0.318 0.3165 0.317 



description of the ^iPIV setup can be found in [9]. The 
concentration of the stained DNA molecules, based on the 
interrogation volume, was less than 8 x 10 particles/ml. 
The images were recorded using a Dantec 80C77 
Hisense PIV 1,344 x 1,024 x 12 bit interface transfer 
camera (Dantec Dynamics A/S, Skovlunde, Denmark). 
A total of five images were taken for each flow field with 
a spatial resolution of 64 x 64 pixels. The interrogation 
cell overlay was 50%. The background noise effect was 
removed by subtracting the background intensity from 
captured images. In addition, an ensemble averaging 20 
images consecutively captured in 4 s was used to obtain 
the velocity measurements and to avoid the Brownian 
motion of the stained DNA molecules. A total of 800 
sets of data were taken at each location for a specified 
Re. The selection of 800 datasets was based on the 
examination of the data convergence. Each measure- 
ment was repeated at least five times under specific 
conditions. 



Results and discussion 

Prior to the formal runs, the velocity in different buffer 
solutions with varied viscosity for the present PZT pump 
should first be calibrated. Through [iPlV measurements, 
average velocity for five different buffers with three dif- 
ferent viscosities of 40, 60, and 80 cP was measured and 
calculated. The results are now plotted against the PZT 
input voltage, as shown in Figure 3. Generally, the distri- 
bution showed a common trend in which a linear pro- 
portionality was present. The higher viscosity caused a 
lower velocity distribution, as expected. The slope of the 
distribution became smaller as the viscosity increased. 
The velocity magnitude spans from 100 to 300 |im/s as 
the input voltage rises from 2.6 to 3.0 V (direct current 
(DC)). The buffer solution effect on the velocity seems 
not to have been noted. 

There are ten semi-circular channels with different 
radii from 500 to 5,000 |im. With different curvature ef- 
fects (i.e., different Dean numbers), the stretching effect 
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Figure 8 Representative of DNA recoiling at different times (At= 5 s) for 1x TBE. 
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differs. It was found that due to the higher Dn, the 
smaller the radius, the longer the stretching. Therefore, 
only data for the radius of 500 \im with Ix Tris-borate- 
EDTA (TBE) and 80 cP at Re = 5 x 10 * (Wi = 12.5) was 
presented, as shown in Figure 4. Seven sequent images 
of the present stretching were illustrated with different 
stretching ratios at the corresponding time. A total 



period of a cycle takes about 9.6 s with each time inter- 
val of 1.6 s. The maximum stretch occurred at the center 
of the semi-circular duct. The stretch ratio was oscilla- 
tory rather than monotonic due to the pressure recovery 
when the flow moved though the curved channels. An 
accompanying plot of the local velocity distribution for 
each stretch was also provided to depict the local 
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Figure 9 Graphs showing (a) relaxation time vs viscosity and (b) \i vs y. 
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velocity gradient. This also explains why the maximum 
stretch occurred at the center of the semi-circle, and the 
oscillating behavior began. As stated previously, the local 
velocity fields developed via i^PIV can be used to quan- 
tify the magnitude of the flow around the semi-circular 
duct, as well as the strength of the shear force. In each 
image, the DNA molecule stretch was clearly observed as 
the corresponding stretch ratio increases, confirming cyc- 
ling between stretched {0<d< 90°) and relaxed (90° <6< 
180°) forms. Due to the parabolic velocity profile, the 
DNA stretch was not uniform across the microchannel 
and DNA molecules near inner walls were more stretched 
than those occupying the central portion and outer wall of 
the channel due to the centrifugal force. 

Figure 5a shows the mean stretch ratio distribution 
versus time in two different buffer solutions with different 
Wi (7.3 to 12.4). As expected, the buffer solution seems to 
exhibit no significant influence on the stretch ratio; it in- 
creases as the Wi increases. In addition, the mean stretch 
seems constant and is independent of time in a time 
period of 6 min. DNA molecule elongation was plotted 
against time and is shown in Figure 5b, in which an expo- 
nential decay form was found for three different viscos- 
ities: 40, 60, and 80 cP. The longest elongation was 
secured with a viscosity of 80 cP, as expected, while the 
shortest is for 40 cP. Taking a close-up look, one may find 
different relaxation times of 3.8, 5.6, and 7.6 s for different 
viscosities of 40, 60, and 80 cP, respectively. With time 
passing, elongation of the DNA molecules reaches a mini- 
mum for each viscosity which has a value of 1.9, 2.2, and 
2.3 |im for the corresponding viscosities of 40, 60, and 
80 cP at a time of about 13 s. 



Figure 6a,b,c depicts the DNA molecule stretch ratio 
histogram for all five different buffers with three viscos- 
ities, respectively, for Wi (Re) from 7.6 (0.3 x 10"^) to 12.5 
(0.5 X 10"^). Generally, buffer dependence again seems not 
to have been noted; furthermore, most DNA molecules 
(about two thirds) are in the range of stretch ratio less 
than 0.2 regardless of the buffers and viscosity, although 
this value (0.2) would increase as the viscosity increases. 
For instance, with the highest viscosity of 80 cP, there were 
about 5% of DNA molecules in which the stretch ratio 
could reach to 0.65. Common features for each among 
these three different viscosities can be seen; it was found 
that the extension was positive, and the minimum stretch 
ratio was approximate 0.1 of 40% to 45% of the DNA 
molecules. The stretch ratio would increase to 0.65 as 
the Wi > 11 for viscosity of 40 and 60 cP, as shown in 
Figure 6a, b; for the viscosity of 80 cP, this happens when 
Wi > 7.6, which can be seen in Figure 6c. In addition, 
more than 5% of the DNA molecules can reach this value 
(i.e., stretch ratio 0.65) for Wi = 12.5. Comparisons were 
also made, as shown in Figure 7, with those related studies 
for the viscosities of 40 and 80 cP. The present data are 
consistently higher than those of previous studies [2,10] 
with regard to both the percentage of the stretched DNA 
molecules and their stretch ratio. In fact, about 10% of 
DNA molecule stretch can reach the ratio of 0.52, and 
about 7% of DNA molecules can reach 0.63. Again, these 
are higher levels than those of previous studies. Table 4 
shows a summary of the DNA mean stretching rate for all 
the cases under study. 

Based on the DNA molecule conformation history, it 
was found that the entire semi-annular duct exhibited 
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Figure 10 Comparisons with those of previous studies. 
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two different opposite trends. First, in the first half duct like recoiling. This is also evidenced by Figure 8, as time 

(i.e., 6 < 90°), the DNA molecules obviously experienced increases with an interval of At = 5 s. Figure 9a,b shows 

stretching; however, for the second half duct (i.e., the relaxation time versus viscosity and the functional 

9O°<0<18O°), it experienced an opposite behavior relationship of viscosity with respectively. Following 
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Figure 1 1 Graphs showing (a) mean stretch ratio vs Wi and (b) mean stretch ratio vs Pe. 
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Figure 9a, one may conclude that the relaxation time 
was a function of y as well. Also included in Figure 9a 
are those from the Rouse/Zimm model and Fang et al. 
[11] for comparison. Good agreement and consistency 
were found. In fact, the present results for the five dif- 
ferent buffers under study were between those of exist- 
ing models. In Figure 9b, the viscosity which was 
correlated in terms of power law with an average power 
of 0.7 was found under different DC voltage inputs. The 
maximum stretch of the stretching force was plotted 
and is shown in Figure 9a with comparisons to those of 
hsted models [12,13]. The data shown strongly indi- 
cated that a small stretching force was needed, as com- 
pared to the existing model with the same stretching 
length. However, the developing trend of the present 
study is the same as those of existing models [12]. The 
viscosity effect for ^ = 40 ~ 80 cP of the present study 
seems not to have been noted as far as the stretching force 
is concerned, as shown in Figure 10. The Freely Jointed 
chin model (FJC) and Wormlike chain model (WLC) can- 
not be compared due to their small values (approximately 
0.12 pN). The stretch length of the present study is clearly 
quadratically dependent on the stretching force. Due to 
the fact that both FJC [12] and WLC [13] models may not 
be applicable, Hsieh and Liu's results [1] were now in- 
cluded in Figure 10 for comparison. It can be seen that 
less stretching force was needed for the present study 
compared to those of Hsieh and Liu [1]. 

Finally, data for mean stretch ratio were correlated in 
a power law from of Wi as x/L^ = 0.17 Wi"^'^^, as indi- 
cated in Figure 11a. Teixeira et al.'s [14] and Smith et al.'s 
[15] results were also included in Figure 11a. Again, the 
present results show a large stretch with a definite Wi. 
Another correlation of mean stretch ratio as a function 
of Pe is shown in Figure lib. A straight line relation was 
found in the form of x/L^ = 5.37 x 10"^ Pe + 0.18, and the 
initial stretch length was obtained as Pe equals zero in 
this study. 

Conclusions 

DNA molecule dynamics in curved (semi-circle, 0° <6 < 
180°) microchannels with different radii for five different 
buffer solutions of Ix Tris-acetate-EDTA (TAE), Ix Tris- 
borate-EDTA (TBE), Ix Tris-EDTA (TE), Ix Tris-phosphate- 
EDTA (TPE), and Ix Tris-buffered saline (TBS) with a 
variety of viscosity such as 40, 60, and 80 cP were exten- 
sively studied for 10 *<Re <10"^ and 5 < Wi <12. The 
major findings drawn are as follows: 

1. Radius effect was significantly noted with maximum 
stretch ratio occurring at the center of the 
semi-circle {6 = 90°) with a radius of 500 |im. 

2. The oscillatory/recovery nature of the present 
stretching behavior was found. 



3. The buffer solution type seems to have no significant 
influence on the stretch ratio, with no viscosity 
effect. 

4. The correlation of x/L^ was developed for 
parameters of Wi and Pe, respectively, with different 
functional relationships. 
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